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Abstract 



An expression for the partial frequency redistribution (PRD) matrix for line scattering in a two- 
term atom, which includes the J-state interference between its fine structure line components is 
derived. The influence of collisions (both elastic and inelastic) and an external magnetic field on 
the scattering process is taken into account. The lower term is assumed to be unpolarized and 
infinitely sharp. The linear Zeeman regime in which the Zeeman splitting is much smaller than 
the fine structure splitting is considered. The inelastic collision rates between the different levels 
are included in our treatment. We account for the depolarization caused by the collisions coupling 
the fine structure states of the upper term, but neglect the polarization transfer between the fine 
structure states. When the fine structure splitting goes to zero, we recover the redistribution matrix 
that represents the scattering on a two-level atom (which exhibits only m-state interference - 
namely the Hanle effect). The way in which the multipolar index of the scattering atom enters into 
the expression for the redistribution matrix through the collisional branching ratios is discussed. 
The properties of the redistribution matrix are explored for a single scattering process for an 
L = — > I — > scattering transition with S = 1/2 (a hypothetical doublet centered at 5000 A and 
500f A). Further, a method for solving the Hanle radiative transfer equation for a two-term atom in 
the presence of collisions, PRD, and J-state interference is developed. The Stokes profiles emerging 
from an isothermal constant property medium are computed. 

Keywords: Atomic processes - line: profiles - magnetic fields - polarization - scattering - Sun: 
atmosphere 

1. Introduction 

The Solar spectrum is linearly polarized due to coherent scattering processes in the Sun's atmo- 
sphere. This linearly polarized spectrum is as rich in spectral structures as the ordinary intensity 
spectrum, but it differs in appearance and information contents (see Stenflo & Keller [1, 2], Stcnflo 
et al. [3]). It is therefore referred to as the "Second Solar Spectrum". A weak magnetic field mod- 
ifies this spectrum through the process of the Hanle effect. This makes the Second Solar Spectrum 
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Figure 1 : Schematic Level diagram of a two-term atom. The lower term is assumed to be infinitely sharp whereas 
the upper term is both radiatively and collisionally broadened. 



sensitive to a field strength regime that is inaccessible to the ordinary Zeeman effect. It therefore 
has the potential to significantly advance our understanding of the Sun's magnetism. In order to in- 
terpret the wealth of information imprinted in the Second Solar Spectrum, it is necessary to develop 
adequate theoretical tools that can later be used for the polarized line formation calculations. 

Stenflo [4, 5] (hereafter S94 and S98 respectively) developed a classical theory for frequency- 
coherent scattering of polarized radiation in the presence of magnetic fields of arbitrary strength 
and orientation. This theory was generalized by Bommicr & Stenflo [6] (hereafter BS99) to include 
the effects of partial frequency redistribution (PRD) in the scattering process. Their formulation 
was restricted to the rest frame of the atom. The transformation to the laboratory frame was 
presented in Sampoorna et al. [7] (hereafter HZ1) for the special case ofaJ=0-s-1^0 
scattering transition. Recently Sampoorna [8] has generalized the classical PRD theory of BS99 
to treat other types of atomic transitions with arbitrary J-quantum numbers. 

In Smitha et al. [9] (hereafter PI) we derived the polarized PRD matrices for a two-term 
atom with an arbitrary L a — >• Lb — > L a scattering transition, taking into account the effects 
of J-state interference between the fine structure components of the split upper term Lb (see 
Figure 1). However, these expressions were limited to the collisionless regime. In the present 
paper we generalize the semi-classical theory of Sampoorna [8] to include J-state interference for 
a two-term atom in the presence of collisions. In Smitha et al. [10] (hereafter P3), a simpler 
version of this theory (see Section 3.2) has already been applied to model the non-magnetic linear 
polarization observations of J-state interference phenomena in the Cri triplet. 

Collisions play a vital role in determining the polarization properties of the scattered radiation. 
For the case of a two-level atom with unpolarized lower level, Omont et al. [11] developed the 
quantum theory of polarized scattering in a non-magnetic medium, including PRD effects. They 
describe in detail the role played by elastic and inelastic collisions. The effects of magnetic fields 
were considered in Omont et al. [12]. An explicit form of the polarized PRD matrix for resonance 
scattering on a two-level atom was derived by Domke & Hubeny [13], based on the work of Omont 
et al. [11], assuming that the lower level is unpolarized. Under the same assumption, a more elegant 
form of the PRD matrix for both the non-magnetic and magnetic cases was derived in pioneering 
papers by Bommier [14, 15, 16] using the master equation theory. The equivalence between the 
QED theory of Bommier [15] and the semi-classical theory was demonstrated in Sampoorna et 
al. [18] (hereafter HZ2) for a J = -> 1 -> scattering transition, and in Sampoorna [8] for an 
arbitrary J a — > J, — >• J a scattering transition. An alternative PRD theory based on the concept of 
metalevels has been developed by Landi Degl'Innocenti et al. [17] for the collisionless case. This 
formulation can also deal with J-state interference in the presence of magnetic fields. 
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In the present paper, starting from the Kramers-Heisenberg formula, we derive the expressions 
for the collisional PRD matrices including the effects of J-state interference for a two-term atom. 
The following assumptions are made: 

1. Infinitely sharp lower term. 

2. Unpolarized lower term. 

3. Weak radiation held limit (i.e., stimulated emission is neglected in comparison with the spon- 
taneous emission). 

4. Hyperhne structure is neglected. 

5. The effects of inelastic collisions that couple the fine structure states are treated approximately 
(see below). 

6. The depolarizing elastic collisions that couple m-states belonging to a given fine structure 
state Jb are taken into account, but are assumed to be independent of the J-quantum number 
for the sake of mathematical simplicity. 

7. Wc restrict our attention to the linear Zceman regime of magnetic field strengths. 

The assumption of an unpolarized lower term is made for the sake of mathematical simplicity, 
but can often be justified when the lower term represents the ground state of the atom. In the stellar 
atmospheric conditions the ground state is generally two orders of magnitude more long lived than 
the excited state, which makes it correspondingly much harder for any ground state polarization to 
survive collisional and magnetic depolarization, as compared with the excited states (see Kerkcni 
& Bommier [19]). We also ignore the induced emission, because in scattering problems it acts as a 
negative absorption and only affects the radiation in the exact forward direction (scattering angle 
exactly zero). The induced emission probability is nearly three orders of magnitude smaller than 
the spontaneous emission probability (see Kerkeni & Bommier [19]). 

The inelastic collisions between the upper and lower terms are treated exactly while the inelastic 
collisions between the upper fine structure states are treated approximately. The inelastic collisions 
between the upper fine structure states (denoted by Tjj b j b ,) manifest themselves in two different 
ways, (i) through a depolarization of state Jb and (ii) through a transfer of alignment and orientation 
between Jb and Jy- 

Since the colliding particles are isotropically distributed around the radiating atom, they destroy 
the alignment and thereby depolarize the levels. Therefore the inelastic collisions that take the atom 
away from the state Jb always depolarize Jb- They also contribute to the inverse lifetime of Jb under 
consideration. In this paper we take into account such inelastic collisions between the fine structure 
states Jb and Jb' by adding these inelastic collision rates (Tjj b j b ,) to the inelastic collision rate 
^i.hJf (where / is the final state). The depolarizing effects of these inelastic collisions are similar 
to the depolarizing effects of elastic collisions. Thus we merge these two effects and define a common 
damping rate 7b for the Jb state. 

The inelastic collisions between polarized fine structure states can lead to a transfer of alignment 
and orientation between them (hereafter referred to as the transfer of polarization). This is similar 
to optical pumping by radiative transitions. The only difference is that the radiative transitions 
between the fine structure states Jb and Jb' are not allowed. Taking account of such transfer rates 
caused by inelastic collisions actually involves formulating the statistical equilibrium equations for 
the concerned states including the atomic polarization of the various states. This is outside the scope 
of our present paper. A formulation of statistical equilibrium equation including these collisions but 
neglecting the redistribution effects in scattering has been presented in Kerkeni [20] and Kerkeni & 
Bommier [19]. They derive the expressions to calculate these rates taking examples of few atomic 
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systems of relevance to the analysis of the second solar spectrum. Our present treatment of inelastic 
collisions is basically heuristic and only takes into account the depolarizing effects of Tjj b j b , . 

The frequency redistribution function that describes the effect of collisions in unpolarized ra- 
diative transfer is the well known type-Ill (or i? m ) function of Hummer [21]. Here we describe 
the matrix generalizations of this standard collisional redistribution function, brought about by 
the magnetic fields and the J-state interference. Using the method described in Appendix C of 
HZ2, we rewrite the PRD matrices in terms of the irreducible spherical tensors for polarimetry. We 
discuss in detail the procedure to identify the multipolar index K, which needs to be assigned to 
the branching ratios that govern the effect of the depolarizing collisions. We illustrate the effects 
of collisions on the Stokes (I,Q/ I,U/I ,V/ 1) profiles of the scattered radiation for the 90° single 
scattering case. Then we present the technique of incorporating this Hanlc redistribution matrix 
for the two-term atom into the polarized radiative transfer equation, and solve it for an isothermal 
constant property atmospheric slab. In the collisionless case, the relevant redistribution matrix 
derived in PI was incorporated into the transfer equation in Smitha et al. [22] (hereafter P2) 
and solved for a constant property isothermal media in the absence of a magnetic field. The same 
method of solution presented in P2 is also used here, but including the collisional redistribution in 
the presence of a magnetic field. 

In Section 2 we derive the elements of the ensemble averaged coherency matrix both in the 
atomic and laboratory frames for a L a — > L b — > L a scattering transition taking into account the 
elastic collisions. In Section 3 we express the type-Ill redistribution matrix in terms of the irre- 
ducible spherical tensors both for the non-magnetic and magnetic cases. The important question of 
identifying the multipolar index K that describes the transfer of angular momentum in a scattering 
event affected by the depolarizing collisions is discussed in detail. The laboratory frame expression 
for the collisional redistribution matrix is also derived in this section. The procedure to incorporate 
this redistribution matrix into the polarized radiative transfer equation for both the magnetic and 
non-magnetic cases is discussed in Section 4. The Stokes profiles resulting from a single 90° scatter- 
ing event and from multiple scattering in an isothermal atmospheric slab are presented in Section 5. 
Concluding remarks are given in Section 6. Finally, in Appendix A we give the expressions for the 
magnetic redistribution functions of type-Ill. 



2. An approximate treatment of collisions including the J-state interference 

2.1. A semi- classical formulation for the polarized PRD matrix 

The Mueller matrix M that describes the transformation from the incident to the scattered 
Stokes vector is given by 

M = TWT \ (1) 

where T and T _1 are purely mathematical transformation matrices. Their explicit forms are given 
in Equation (9) of S98. The W-matrix is defined in Equation (7) of PI. The elements of this 
matrix contain the bilinear products of the complex probability amplitude w a 0(Jffj,fJ a fi a ). These 
amplitudes for transition from an initial state a to final state / via all intermediate states b are 
given by the Kramers-Heisenberg formula (see S98 and Sampoorna [8] for historical accounts) as 

V- {/jf ; gqJ6) (b\r . e p \a) 

w a p ~ > : jz , (2) 

*-f LU b f - W- 176/2 
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where uj = 2ir£ is the angular frequency of the scattered radiation in the atomic rest frame, huib / 
is the energy difference between the excited and final states, and 7b is the damping constant that 
accounts for the broadening of the excited state b, while the initial and the final states are assumed 
to be infinitely sharp. The damping parameter is assumed to be the same for all the magnetic 
substatcs of the excited state. The matrix elements appearing in Equation (2) can be expanded 
using the Wigner-Eckart theorem as 



w a0 (J f fi f J a [i a ) ~ ]T (-i) q - q '^(2J a + l)(2J / + l)(2J b + 1)(2L„ + 1) | L f b L j b f \ \ 

where fib represents the magnetic substates of the upper state b with total angular momentum 
quantum number Jb, orbital angular momentum quantum number Lb, and spin S. The quantities 
J a and Jf are respectively the total angular momentum quantum numbers of the initial and final 
states a and / with orbital angular momentum quantum number L a , and magnetic substates [i a 
and fif. The quantities e are the geometrical factors (see Equations (2) and (27) of S98), with a 
and (5 denoting the outgoing and incoming radiation, respectively. In Equation (3), q — fif — fib and 
q' = fjb a — fib- In the rest of the paper we denote the indices as follows for the sake of convenience 

Jb =b, J a = a, Jf = /; JbfJ-b = b m , JafJ-a = a m , JffJ-f — fm, 

Jb'fJ-b' = b' m , Jbflb" = C> J b'^b"' = C- ( 4 ) 

The frequency-normalized profile function is given by 

®i b {vb m f m - = 1 (^ 1 L 77IZ\ with Vh ^S m = "bf + - 9flif)"L- (5) 

Here hvbf is the energy difference between the upper (Jb) and lower (Jf) states in the absence 
of magnetic fields, gb, 9} are the Lande factors of these states, and vl is the Larmor frequency. 
Equation (3) refers to the case of frequency-coherent scattering in the atomic rest frame. 

The phenomenological extension of Equation (3) to the case of PRD is achieved by treating each 
radiative emission transition between magnetic substates fib and [if in terms of a damped oscillator 
that is truncated by collisions (see HZ1). In other words, in Equation (3) we make the following 
replacement for the profile function : 

*TkK»/m - — > (h m )a mfml (6) 

where the Fourier-transformed solution of the time-dependent oscillator equation is given by (see 
BS99) 

(h m )a mfm = (0»m/ m + C ^bDa mfm . (7) 

Here we have omitted the unimportant phase factor, as it vanishes in the bilinear product 
{^b m )a m f m (fb' m )a m f m - The constant C in Equation (7) defines the relative amplitudes of the sta- 
tionary and the transitory parts of the solution, which are given by 

{K:% m f m = ^b mam - em VamJm ), w 
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)a m f„ 



= %M mam a^M m fm -oh- e- i( ^-^ 2 -^] . (9) 



Here £' denotes the frequency of the incoming photon in the atomic rest frame, t c is the time 
between two successive collisions, co bm f m — ^ v b m f m - The profile function & lb (v bm a m — £') 1S given 
by Equation (5) with £ replaced by while v bm f m is replaced by v bm a m that is defined similar 
to Equation (5). In Equation (8), v am f m appearing in the delta function is demanded by energy 
conservation (see Equation (9.10) of S94), and is given by 

v a m f m = Vaf + (9aHa ~ 9ffJ-f)^L, (10) 

where hu a f is the energy difference between the states J a and J/ in the absence of a magnetic held. 

2.2. Coherency matrix in the atomic rest frame 

The elements of the ensemble averaged coherency matrix (f bm f^, ) 0m / m can be derived starting 
from Equations (8) and (9), applying the same steps that are described in detail in BS99. These 
elements are contained in the bilinear product w a p(f m a m )w* a , p,{f m a m ). In the atomic rest frame 
ensemble averaged coherency matrix elements are given by 

(r bm r* b >Ja mfm - A W cos^ m6m e^'^ *Z^f- m (0<*(£ ~ ? ~ «WJ 

+B W co S (3 b , mbm cosa b , mbm e^^+^^^ (11) 

where the angles flb' m b m and a b i mbm (arising due to the combined effects of the J-state and m-state 
interferences) are defined respectively by 



, a ^vb + {gvw - 9bHb)u L , u Vb + (g b >fi b < - g b fi b )u L 

t&nfj b , mbm = ■ ; tanav m 6 m = ■ — , 

ivb + 7c lb' b + 7c/2 

with -f bb > given by 



(12) 



lb + W n 
lbb> = ^ = 7b ' b ' ( 13 ) 

Here j c is the collisional damping constant, while fuv b > b is the energy difference between the J b > and 
J b states in the absence of a magnetic field. The elastic collisional rates are in general different for 
each fine structure component ( J b ) of the upper term. However, for simplicity we assume them to 
be independent of the J-quantum numbers. 

A bb i and B bb / are the branching ratios for a two-term atom. The explicit expressions for them 
will be defined later in Section 3.1. 

Like in PI, we limit the treatment to the linear Zeeman regime, in which the Zeeman splitting 
is much smaller than the fine structure splitting. When J b ^ J b i the contributions from the second 
terms with cjl in Equation (12) to the angles (3 b > mbm and a b > mbm can therefore be ignored, because 
they are insignificant in comparison with the first terms. The classical generalized profile function 
is defined as 

*£& m / m (0 = \ K»A» - + K> (^m/m ~ 0] - 

in the same way as in BS99. 
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2. 3. Coherency matrix in the laboratory frame for type-Ill redistribution 

We transform Equation (11) to the laboratory frame using the same steps as described in 
Section 2.2 of HZ2. Thus the ensemble averaged coherency matrix in the laboratory frame is given 

by 



(h m r*b>Ja mfm = A w cos)8 6 , m6m e i ft'm*m \{ti b l mfi , m ) amfm + i(/£, i6 ,J„ m /, 



+Bbb> cos (3b' m b m cos ab' m b m e 



,111 , • 

n b m a m ,b' m f m +l h m a m .b' m f, 



III 



(15) 



The various auxiliary quantities for type-II redistribution are defined in Section 3 of PI. Hence 
we do not repeat them here. Hereafter we confine our attention to the collisional redistribution 
(type-Ill). The corresponding derivation for pure radiative (collisionless) redistribution (type-II) 
are given in PI. The auxiliary quantities for type-Ill redistribution that appear in Equation (15) 
are defined by 



h m 

n b m a m ,b' m U 



pill, HH R III,HH pIII.HH , r>i 

n b' m a m .b'mf m + U b' m a m ,b m f m "+" U b m a m ,b'mfm + U b m a m ,b m f„ 



-,111, HH 



R III,FH pIII.FH pIII.FH 

U b' m a m ,b' m f m + U b' m a m ,b m f m U b m a m ,b' m f„ 



R 



III,FH 

fr?n. Q>m fom fir 



■ (16) 



Similarly we have 



jrlll 

lb m a m ,b' m f„ 



1 

4 

+- 



pIII.HF _ dIII,HF i pIII,HF t-,111. HF 

fm b' m a m ,b m f m ~ b m a m ,b' m f m b m 

O'ln i^m frr 



-.III, FF 



nlll.FF _ pill, FF _ pin, FF p i 

U b' m a m ,b'mfm n b' m a m ,b m f m U b m a m .b' m f m + n b m a m .b m f„ 



(17) 



The magnetic redistribution functions of type-Ill appearing in the above equations are defined in 
Appendix A. 



3. The PRD matrix expressed in terms of irreducible tensors 

The importance of expressing the PRD matrices in terms of the irreducible spherical tensors 
introduced by Landi Degl'Innocenti [30] has been discussed in PI. The definition and properties of 
irreducible spherical tensors are described in detail in Landi Degl'Innocenti & Landolfi [23] (here- 
after LL04). The way to incorporate these tensors in the analytic form of the PRD matrix derived 
from a semi-classical approach has been described in HZ2 (see also Section 4 of PI). Applying 
the same method we have obtained an expression for the type-Ill redistribution matrix in terms of 
irreducible spherical tensors. The case of the type-II redistribution matrix has been discussed in 
Section 4 of PI. 

As in PI we now express the type-Ill PRD matrix derived in Section 2 in terms of Tq{i : n), 
where i = 0, 1, 2, 3, and K = 0, 1, 2 with — K < Q < +K. Following the same procedure as discussed 



7 



in Section 4 of PI, the matrix Tj* of Equation (22) in PI, which describes the transformation of 
the elements of the coherency matrix, can be written in the atomic rest frame as 

T^ pa (Z,n;Z',n',B) = (2L a + I) 2 £ G Z 6 Z 3 (-l)«-«'+«"-«"'f|,„( / i,^n)^, v (a,p,n') 

a m fmb m b' m 

x{A bb , co S p b , mbm e**^ *£££ m (0<S« " ^ " "« m /J 

+B W cos/3 b , mbm cosa^e^^+^Oc^^^^ (18) 
where 

G = (2J a + 1)(2J / + 1)(2J 6 + l)(2J b , + 1), (19) 

*-{ j J' ill 5" j 1 ill J* 5' JUS* j il- < 20 » 

[ Jb J/ 5 J [ Jb Ja O J [ Jb> Jf S J ( Jb' * J 

and 

( Jb Ja i \( Jb j f i \( j v Ja i \( j v jf i #/ y (21) 

In Equation (18), £| „(/i, i/, n) is a reducible spherical tensor. After transforming to the Stokes 
formalism (see Section 4 of PI), the redistribution matrix for J-state interference can be written in 
symbolic form as 

Ry(£, n; n', B) = n; n', B) + R™(£, n; ?,n',B), (22) 

where the pure radiative part of the redistribution matrix is given by branching ratio A bb > times 
Equation (25) of PI, and the collisional frequency redistribution is taken into account through 

B$ I {€,n;?,n',B) = l(2L a + l) 2 ]T G Z 6 B bb ,^(2K> + 1)(2K" + 1) 

K'K"QafW 
1 1 K" \ ( 1 1 K' 

q -q" Q ) { q'" -q' -Q 



x| (-v q " +q,+Q 

X \ [*7»+7c("hma m - O + $ 76 ,+ 7 >6' m a„, - O] [*7 6 +7>W' m - + $ 7 „, + lc (Vf m f m ~ 

xcos/3 b , m b m cosa 6 , m6 „, e i ^»'m»m+ 01 »'m*m)|(-i)07^"'(i,n)7^(j,n / ). (23) 

Note that in the formal expression for R the branching ratios are built into the R n and R m 
components. As the collisional branching ratio B bb > depends on index K, our next task is to 
determine the explicit form of this dependence. This will be done in the next subsection. 

3.1. Identification and physical significance of the multipolar index K in the collisional branching 
ratios 

It is well known that the spherical unit vectors form a natural basis to decouple the classical 
oscillator equation. Fano [24] suggested that a convenient basis to be used when dealing with 
scattering problems in quantum mechanics, is the irreducible tensorial basis instead of the standard 
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I JM) basis of Hilbert space. This is due to the fact that irreducible tensors transform under co- 
ordinate rotations like the spherical harmonics (Yj m ) and are thus suited for a study of rotationally 
invariant processes. With irreducible tensorial operators one can express the scattering matrix such 
that it formally looks the same in the magnetic (with the polar z-axis along B) and the atmospheric 
(with the polar z-axis along the atmospheric normal) reference frames. This is the advantage of 
going to the irreducible tensorial basis (hereafter called the KQ basis). A more detailed historical 
background for the irreducible tensorial operators is given in Sahal-Brechot et al. [25]. 

Thus the geometrical factors associated with the scattering problem, and also the density matrix 
for the atomic levels in question, should be transformed to the KQ basis. The transformation of 
the geometrical factors to the KQ basis is described in Chapter 5 of LL04 and is used in HZ2. 
The density matrix is first written in the standard | JM) basis and then transformed to the KQ 
basis (see Equation (3.97) of LL04), which is then called 'multipole moments' of the density matrix, 
or 'irreducible statistical tensors'. In the case of the radiation field, the multipole index K has 
the following interpretation : K = means isotropic scattering, K = 1 is related to the circular 
polarization, while K = 2 is related to the linear polarization. In the case of the atomic levels, 
K = represents the population of the level under consideration, K — 1 is related to the orientation 
of the atom, while K = 2 is related to the alignment of the atom (this physical interpretation can 
be found in pp. 128 and 129 of LL04, Section 10.4 of S94, and inTrujillo Bueno [26]) 

In the case of radiation field an irreducible tensor Tq(i, n) is constructed by forming a suitable 
linear combination of the direct product of two geometrical factors. Since geometrical factors 
basically contain unit polarization vectors of rank one, their direct product represents a second 
rank tensor, with K taking values 0, 1, and 2. Note that these values of K can also be obtained 
through angular momentum addition of two tensors of rank 1. In the case of the density matrix of 
the atom, the value of K is determined by the addition of angular momenta J and J'. For example 
for a two-level atom with unpolarized ground level the value of K relating to the statistical tensor 
of the upper level is given by angular momentum addition of J\> and Jw ■ Further, Q takes values 
—K to +K in steps of one, and is related to the magnetic quantum numbers of the upper level. 

We can denote K' as the multipole component of the incident radiation field, K as the multi- 
pole moment of the upper level of the atom, and K" as the multipole component of the scattered 
radiation. The scattering process can be understood as a transfer of the K' multipole component of 
the incident radiation to the K multipole moment of the atom's upper level through an absorption 
process, followed by a transfer of the K multipole moment of the atom's upper level to the K" 
multipole component of the scattered radiation through spontaneous emission. The depolarizing 
collisions that govern the branching ratios and the magnetic field that governs the Hanle angles 
affect the upper level of the atom directly and modify the K multipole moment of the atom, but 
they influence the scattered radiation only indirectly, through spontaneous emission from the level 
that has been directly affected. Thus it is the K index of the upper level of the atom that needs 
to be assigned to the branching ratios and the Hanle angles, and not the multipole component of 
the incident or the scattered radiation. In the absence of magnetic fields or in the presence of weak 
magnetic fields (Hanle effect), K' = K = K" . In the presence of a magnetic field of arbitrary 
strength (Hanlc-Zeeman regime) all the three if's are distinct. This is due to the distinction pre- 
served through the profile functions, which become different for the different Zeeman components. 
However, in weakly magnetic cases (when the Zeeman splitting is much smaller than the effective 
line width) the distinction is so small that it can be ignored. 

From the above discussion it is clear that for a correct identification of K for the branching ratio 
we need to know the density matrix of the upper level in the KQ basis. Since the density matrix 
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does not appear directly in the Kramers-Heisenberg approach that we use, we need to indirectly 
identify K either by drawing analogy with the density matrix theory, or by using a suitably defined 
quantum generalized profile function. Such a function was defined by Landi Degl'Innocenti et al. 
[27] for the special case of a two-level atom (without J-state interference). The multipolc moment 
K of the upper level is built into this function through the third 3 — j symbol appearing in the 
following definition : 

*% K '(Ja, J b ; O = V3(2J6 + 1)(2K + l)(2K> + 1) {-l) Jb ~^- 1+Q 



fJ-bfJ-b' fJ-aPP' 



Jb Ja 1 

~W> Ma P 



Jb Ja 1 
~Mb' Ma p' 



Jb K J b 

-Hb Q Mb' 



1 1 K 1 

-P P' Q 



(24) 



The $q' K defined above can be seen as a frequency-dependent coupling coefficient that connects 
the (K , Q) multipole component of the incident radiation field with the (K, Q) multipole moment 
of the atomic density matrix (see LL04, p. 525). In the non- magnetic and weak field limits, the vl 
dependence of the profile function <& 7() + 7c can be neglected, which gives us 



(25) 



where is defined in Equation (10.11) of LL04, and <f> denotes the usual non-magnetic profile 

function. In this limit we have K' = K. 

In the case of J-state interference a suitable quantum generalized profile function has not been 
defined yet, but we can define it here in analogy with the two-level atom case. It has the following 
form for the incoming radiation: 

<$>Q K '(a, b', b; O = (2 J a + l)y/3(2J b , + 1)(2J„ + l)(2K + 1)(2K< + 1) ]T (-l)^"^' 

MbMb' Vaq'q'" 

Jb Ja 1 \( Jb' Ja 1 \( Jb Jb' K \f 1 1 K' 

-fi b fi a -q' J \ -fi b ' Mo -q'" / v Mb -m^ -Q ) v i'" -i' ~Q 

1 



COS (3b' m b m £ 



(26) 



with a similar expression for the outgoing radiation when J a and fx a are replaced respectively by J/ 
and Mf, an d angle fib> m b m is replaced by oib' m b m - Notice that unlike the two-level atom case we now 
have included the angles /?b' m b m an d ab' m b m in the definition of the quantum generalized profile 
function, as they cannot be taken outside the summation over the magnetic substates. Using the 
orthogonality relation of the 3 — j symbols, it is easy to verify that 



J2 ®q K («, b', b- 0$q' K (/, b', b- = 3G y/(2K' + l)(2K" + 1) 



K 

X 



1 

X - 

2 



E 



Zz (-l)9'+«"+0 

q'" 

^1 b +-iA l/ b m a m - £,') + $ 7b , +7c (^b' m a„ 
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cosP b , mbm coaa b , mbm e li - li -'^ +a -'^\ (27) 



which is a useful relation that helps in the identification of the multipolar index K . An equivalent 
relation, but for the case of m-state interference, is Equation (22) of Bommier [15]. 

Comparing the terms in the flower brackets of Equation (23) with the RHS of Equation (27), 
we can see that they are the same (except for some factors). Therefore after substituting the terms 
in the flower brackets of Equation (23) with the LHS of Equation (27), we assume 7 c /2 = as 
a reasonable approximation (see S94), where D^ K "> is the 2K multipole collisional destruction rate. 
Further, following BS99, we identify 7& = Tj^ b + Ti b and 7 C = Te, where r Rb is the radiative width 
of the fine structure state J b . Tj b is the total inelastic collision rate defined for the state J b . It is 
given by 

r Ib = Y, T ibf + Y, Tlbb '- ( 28 ) 

Here Ti b f couple the upper state J b to the lower state J/ and Tj bb ' couples the two fine structure 
states J b and J b > . Indeed such a definition of total inelastic collision rates can be found in Omont et 
al. [11] and also in Equations (2.15)-(2.20) of Heinzel & Hubeny [28]. Te is the elastic collision rate 
and D^ K } represent the depolarizing elastic collisions that couple the Zeeman substates (m-states) 
of a given ^-statc. In general D^ K "> may be different for each of the fine structure components 
with quantum number J b . However, as an approximation we assume them to be independent of 
the J-quantum number. 

With the above mentioned substitutions and identifications, Equation (23) can be rewritten as 

n; n', B) = 2 -(2L a + l) 2 £ Z 6 {-\)*Tf (i, n)T K Q (j, n') 

K K' K" Qafbb' 

x^ K '(a,b',b;^ K "(f,b',b^), (29) 
where B^f? is the collisional branching ratio defined as 

Also, the branching ratio A bb > can be written as 

A »»' = w Jb' r . ( 31 ) 

L R +1/ +L E 

where 

■pbb' _ F Rb + T Rb , —66' _ Tib + Ti b , 

l r — 2 ' 1 — 2 ' 

with T/6 defined in Equation (28). The total damping rates that appear in the branching ratios 
are the same as those that appear in the denominator of the Hanle angles (see Equations (12) and 
(36)). Therefore the J b J b > dependence of the branching ratios defined now for a two-term atom is 
self-consistent. We have verified that when we set J b — J b > and J a — Jf (the case of a two-level 
atom with only m-state interference) in Equation (29), we recover Equation (49) of Bommier [15]. 
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3.2. The laboratory frame expression for the redistribution matrix R m 

We convert Equation (29) into the laboratory frame using the same procedure as described in 
Section 2.2 of HZ2. The resulting expression for the normalized type-Ill redistribution matrix in 
the laboratory frame can be written as 



R]f(x,n;x',n',B) = 



2L b + l x ^ 

25*+ 1 ^ 

KK'K"Qafbb' 



Z 6 Bi«\-l)QTf {i,n)T« Q {j,n') 



(33) 



where T^q\^' K {x, x' , O, B) is the laboratory frame redistribution function obtained after transfor- 
mation of the atomic frame functions $g' K '(a, &', b; £')$q' K " (f, b' , b; £)• The factors (2L b +l)/(2S+ 
1) result from the renormalization of Equation (29). The function TZq y^' K has the following form: 



Kqui {x,x,Q,B) = 



x J h -W„+q CQS/ 
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(34) 



The results presented in Section 5 are computed using Equation (33) of the present paper for the 
type-Ill redistribution and Equation (25) of PI for the type-II redistribution matrix. The form of 
the total redistribution matrix in the laboratory frame remains the same as Equation (22). In the 
non-magnetic case Equation (33) simplifies to 



R 



in 



(x,n;x',n')= 3 -^^- £ G Z 6 

(-!)%*(;, n)T5 Q (j,n'). (35) 



Lin , ■ fin 

n ba,b'f + l Jba,b'f 



KQafbb' 
1 1 K 
Jb' Jb Ja 



1 1 K 

Jb' Jb Jf 



The angle fob an d the auxiliary functions h b ^ b ,j and fb a \'f are defined respectively in Equa- 
tions (12), (16) and (17), but with v L = 0. The angles and fob (arising exclusively from 
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J-state interference) are denned as 

, (K) ^b'b , a ^b'b /oa , 

tana),,' = —rrr ttt ; tanrVb = —rrr tjt , (oo) 

with T R — and b = (see Equation (32)). In P3, the two-level atom branching ratios 
were used in the realistic modeling of the linear polarization profiles of the Cri triplet. These 
branching ratios can be recovered from the more general two-term atom expressions given in Equa- 
tions (30) and (31) by neglecting Tibb'- This is equivalent to setting Jb = Jv in Equations (30) and 
(31). The angle- averaged redistribution matrices corresponding to the angle-dependent redistribu- 
tion matrices presented in Equations (33)-(35) can be recovered by replacing the angle-dependent 
redistribution functions (Equations (A.1)-(A.4)) by their angle-averaged analogues. These angle- 
averaged functions are obtained by numerical integration of the angle-dependent functions over the 
scattering angle (see Equation (A. 6)). 



4. The polarized radiative transfer equation for J-state interference 

The polarized radiative transfer equation for the Stokes vector J in a one-dimensional planar 
medium for the Hanle scattering problem can be written as 

/ I(T ^' n) = Wx)+r)[I(r,x,n) - S(r,x,n)}, (37) 

where the notations are the same as those used in P2, with the positive Stokes Q representing electric 
vector vibrations perpendicular to the solar limb. This definition is opposite to the way in which 
the positive Stokes Q is defined in the observed spectra. This can easily be accounted for (through 
a sign change), when comparing the observed spectra with the theoretical results, n = (~d,(p) 
defines the ray direction where $ and ip are the inclination and azimuth of the scattered ray with 
fi = cos# (see Figure 1 of PI). In the weak magnetic field limit, the Stokes vector I = (I, Q, U) T 
and the Stokes source vector S = (Si, Sq, Sjj) T ■ In this limit, the transfer equation for Stokes V 
decouples from that of the Stokes vector (I, Q, U) T . This is known as the weak field approximation. 
In Equation (37), the Stokes vector I and the Stokes source vector S depend on n. In the case of 
angle-averaged redistribution, it was shown by Frisch [29] that one can decompose S and I into six 
cylindrically symmetric components Iq and Sq with the help of the irreducible spherical tensors 
for polarimetry (See Landi Degl'Innocenti [30]). Here, K = 0, 2 and —K < Q < +K. Such a 
decomposition results in a reduced Stokes vector X which is independent of <p and a reduced source 
vector S which is independent of both i9 and ip. We denote the quantities in the reduced basis by 
calligraphic letters and in Stokes basis by Roman. In such a reduced basis the transfer equation 
can be written as 

mj^v) = {4){x) + r) [x(T; ^ ^ _ s ^ x)] (3g) 

The reduced source vector is defined as 

S{T, X )= * {x)S f*\ + rG{T \ (39) 
4>{x) + r 

where Q[r) = {B, 0, 0, 0, 0, 0} T is the primary source vector. The reduced line source vector is given 

by 



Si{r,x)=Y,VwQ{r)+ I °° nbA *; X '' B) J{T,x>)dx> 



(40) 
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where Hbb' {x, x' , B) is the redistribution matrix for a two-term atom, with the summation over 
Jb and Jf not yet performed. The thermalization parameter is given by 

rf 

£ bb' = by 1 bb > ■ ( 41 ) 
^ R + * / 

The computation of the above defined reduced line source vector is very expensive because of the 
summations over Jb and Jf,/ which need to be performed at each iteration. However for all practical 
applications, we can assume e to be the same for all the JbJv states, which is a good approximation. 
Such an approximate e is constructed by taking an average value of Tjb for all transitions involving 
Jb, Jb' and Jf and an average value of Tr for all the upper fine structure states. Under such an 
approximation, the reduced line source vector can be written as 

r+°° iz( T r > b) 

Si{t,x) = eg{r) + / ^±l£l j(r,x')dx'. (42) 

J-oo <p(x) 

The mean intensity J{r, x) is defined by 

J(r,x) = \j + ^ *(//)Z(T,z,//)d//. (43) 

The elements of the *f?([i) matrix are given in LL04 (see also Appendix A of Frisch [29]). lZ(x, x' , B) 
appearing in Equation (40) is a (6 x 6) diagonal matrix. The explicit form of this redistribution 
matrix with and without the presence of magnetic fields is defined in the following sections. In the 
absence of a magnetic field, only the 2q and Iq components contribute to the Stokes vector. Hence 
the (6 x 6) problem reduces to a (2 x 2) problem. The transfer equation defined in Equation (38) 
is solved using the traditional polarized accelerated lambda iteration technique presented in P2. 

4-.1. The redistribution matrix 7L(x,x') for the non-magnetic case 

In the absence of a magnetic field the redistribution matrix in Equation (40) becomes indepen- 
dent of B and reduces to a (2 x 2) diagonal matrix with elements TZ(x, x')=diag (71°, 1Z 2 ). The 
elements 1Z K are defined as 

k k m = ^\ 1) v^6(-i)W ] ) M/ 1 ) f } 

V ' ' 25+1 ^ K ' \ Jb' Jb Ja J 1 Jb' Jb Jf J 

bb'af 

X[A W COS(3 b 'b ^ b ' b [{hlb')af+Vl\')af 

+<> cosh'b cosagf e^-^) [h™ b , f + if™,,] }. (44) 

The {h\ l b ,) a f and (/ 6 n b /)a/ are the auxiliary functions for type-II defined in Equations (14) and (15) 
of PI and the auxiliary functions for type-Ill are defined in Equations (16) and (17) of the present 
paper. They are used here for the non-magnetic case and with the angle-averaged redistribution 
functions of type-II and type-Ill. In the limit of a two-level atom model (Jb = J' b and J a = Jf), 
the {h l b l b ,)af and h l b l b ,, go respectively to R u and R m functions of Hummer, whereas the (f b u b ,) a f 

and f™v f an< ^ tne angles j3 b 'b and a^ b go to zero. 
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4-2. The redistribution matrix TZ(x,x',B) for the magnetic case 

The redistribution matrix for a two-term atom denned in Equation (22) involves summations 
over the total angular momentum quantum numbers and the corresponding magnetic quantum 
numbers. This does not allow direct decomposition to go from the Stokes vector basis to the 
reduced basis (see Section 4 for details on these two basis) . Such a decomposition is possible in the 
non-magnetic case. This is because in the absence of a magnetic field, the summations over the 
magnetic quantum numbers can be analytically performed using Racah algebra as shown in PI for 
type-II redistribution and Equation (35) for type-Ill redistribution. However in the magnetic case, 
all the summations remain intact and have to be performed numerically. This is very expensive. 
Because of these difficulties, we need to resort to the weak field approximation which allows us to 
apply the decomposition technique. In this regard, the summations over the J-quantum numbers 
can be split into three different terms namely 

n(x,x',B)=J2 K^ a (x,x',B)+Y J nl b ,^ f (x,x',B)+Y J ^b,a,f(x,x',B). (45) 

b=b',a=f b^b',a,f b=b'M^f 

The first term represents the case of 'Resonance' scattering in a two-level atom model with a 
summation over all the lines of the multiplet (see Figure 2a). This contributes mainly to the cores 
and near wings of the lines within the multiplet. Its weak field analogue has already been derived 
in Bommier [15] and can be expressed as 

K A (x, x', B) = Ka& B ) = E Wb,aK.a(x, *' , B). (46) 

ab ab 

Here 7t^ a (x, x' , B) is the Hanle redistribution matrix for a two-level atom with J a — > J — > J a 
scattering transition as presented in Bommier [15]. This is also the same redistribution matrix 
defined in Appendix A of Anusha ct al. [31], but for a J a — > J — > J a scattering transition. The 
details of the domain based decomposition of this matrix are also given in the above paper. In the 
reduced basis, 1Z A (x, x' , B) is a (6 x 6) matrix. Wb, a are the weights for each line component of 
the multiplet (derived from Equation (44) with Jj, = Jy and J a = Jf) and are given by 

w »- - 1 irr< 2J " + 1 ) 3 < 2J ' + 1 > 2 {!;; i 1 1} 1 - < i7 > 

The second term represents only the J-state interference between different lines of the multiplet 
(see Figure 2b). It includes both the 'Resonance' and the 'Raman' scattering parts and is effective 
mainly in the wings between the lines. This term is quite insensitive to the strength of the magnetic 
field. This can be seen from Figures 2 and 3 of PI where the magnetic field effects are confined 
mainly to the line cores. Hence in this component we can set the magnetic field equal to zero as a 
good approximation. This makes the evaluation of both the Resonance and the Raman scattering 
parts similar to that performed in P2 and P3. In the reduced basis this term is simply given by a 
(6 x 6) matrix lt B (x,x') which is equal to diag(7\L°, 1Z 2 , 0, 0, 0, 0). Here 1Z K are the redistribution 
functions which include the effects of collisions and the J-state interference between different line 
components in a multiplet defined in Equation (44). Only J ^ J' b terms are retained in the 
summations appearing in this equation. The J b = Jv contributions are contained in the first term. 
However in some of the well known examples in the second solar spectrum like the Mg n h and k, 
Call H and K and the Cri triplet, the initial and the final states are the same. Also for the case 
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Figure 2: The schematic level diagrams representing the three components of the J-state redistribution matrix. 
Panels (a), (b), and (c) represent respectively the two-level atom resonance scattering, J-state interference, and 
Raman scattering (resonance fluorescence). 



of the hypothetical doublet considered in this paper, arising due toanL = 0->1^0 scattering 
transition with spin S = 1/2, the initial and the final states are the same. Hence the Raman 
scattering part does not play a role. 

The third term represents the case of only Raman scattering without the J-state interference, 
where the initial and the final states are different (see Figure 2c). A derivation of the weak field 
analogue of this component (in a way similar to that of Bommier [15]) is yet to be performed. 
Again for some of the well known examples mentioned above, this component does not contribute. 
Thus the final expression for the redistribution matrix that is used in Equation (40) is 

TZ(x, x', B) « Tl A (x, x', B) + K B (x, x'). (48) 



5. Results and discussion 

In this section, we study the effects of collisional redistribution matrix on the emergent Stokes 
profiles for the case of single scattering and also multiple scattering in an isothermal atmospheric 
slab. All the profiles presented in this paper are computed for a hypothetical doublet line system 
with the line center wavelengths at 5000 A and 5001 A arising due to an L = - > 1 — > scattering 
transition with spin S = 1/2. The J quantum numbers of the lower and upper states are J a = Jf = 
1/2 and Jb = 1/2,3/2. In Section 5.1 we present the scattered Stokes profiles resulting in a single 
90° scattering case. In Section 5.2 we present the multiply scattered Stokes profiles emerging from 
an isothermal constant property atmospheric slab with and without the presence of a magnetic 
field. 

5.1. The single 90° scattering case 

To explore the general behavior of the redistribution matrix in the presence of collisions we 
illustrate the Stokes profiles that result from single 90° scattering event. We examine the influence 
of the elastic collisions on the Stokes profiles in the presence of a magnetic field. The magnetic field 
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Figure 3: Effect of collisions: The profiles of the intensity / and the fractional polarizations —Q/I, —U/I and V/I are 
plotted for a hypothetical doublet at 5000 A and 5001 A with field strength parameter «jj = 0.004. The coherence 
fractions used are ^/coh — 0.9 (thick solid line), 'jcoh — 0.5 (dotted line), and *ycoh — 0.1 (dashed line). The fine 
structure splitting is 1 A. Single 90° scattering is assumed at the extreme solar limb (fj, = 0). The model parameters 
are a R = 0.001, i? s = 90° and <p B = 45°. The Doppler width AA D = 0.025 A. 



orientation is given by -&b = 90° and ips = 45° where the colatitude $b and azimuth ipB characterize 
the magnetic field orientation with respect to the polar z-axis (see Figure 1 of PI). We consider an 
unpolarized (I in = [1,0, 0,0] T ) and spectrally flat (frequency independent) radiation field that is 
incident in the vertical direction (parallel to the polar z-axis). The singly scattered Stokes vectors 
are then exclusively determined by the first column of the angle-dependent redistribution matrix 
by integrating over the incident wavelengths. However in the multiple scattered solutions discussed 
in Section 5.2, we restrict our attention only to the angle-averaged redistribution matrix. The 
magnetic field strength is parametrized by the splitting parameter vn given by 

where B is the magnetic field strength, eo is the charge of the electron and m its mass. AAd is the 
Doppler width and is assumed to be 0.025 A for both the lines. The radiative width of the upper 
state is parametrized as am, = T jy, / ' (47T Az^d) • It is assumed to be the same for both the lines and 
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is chosen to be 0.001. The radiative width anb is related to the total damping parameter through 



We assume the inelastic collision rate Tn, to be zero. 

The depolarizing collisional rates = 0.5Te, and = 0. For simplicity we set = 
Z?( 2 '. However in general they can differ (for example, £>W = 0.43r E and = 0.38r E according 
to Berman & Lamb [32]). We have verified that the Stokes V/I is insensitive to the values of 

The collisional effects are built into the R matrix derived in Section 3 through the branching 
ratios defined in Equations (30) and (31). The elastic collision rate is parametrized through the 

coherence fraction 7coh as Icoh 

= l/[l + {T E /T R )}. In the present paper we take Tr to be the same 
for both the upper fine structure states and j C oh to be the same for all JbJf combinations. When 
"icoh — 1 the redistribution is entirely radiative (only R 11 ), whereas jcoh — represents purely 
collisional redistribution (only R m ). Wc consider a range of values jcoh G [1,0] to represent an 
arbitrary mix of R 11 and R m type redistribution. 

Figure 3 shows the Stokes (I, Q/I, U/I, V/I) spectra for a doublet. The polarization of the line 
at 5001 A is zero because its polarizability factor Wi = 0. The collisions affect the wavelength 
domain outside the line core region of this line. But for the line at 5000 A, the collisional effects arc 
seen both in the line wings and the line core. In the line core the collisional effects compete with 
the Hanle effect and in the wings it is an interplay between the J-state interference effect and the 
collisional redistribution effect. 

The value 7 co /, = 0.9 corresponds to a mix with 90% of R 11 and 10% of R m (see the thick solid 
line in Figure 3). The profiles look similar to those for pure R n (see dotted line in Figure 3 of 
PI). However, in Q/I there is a small depolarization, mainly in the wings, due to the presence of 
collisions. The core of the line at 5000 A seems to be less affected than its wings. The Q/I at the 
5001 A line remains zero. In the presence of elastic collisions, a small U/I signal is generated in the 
wings of the 5000 A line. This non-zero U/I wing polarization and the depolarization in the wings of 
Q/I are induced by the elastic collisions in combination with the magnetic field and can together be 
referred to as the 'wing Hanle effect'. This effect arises because the elastic collisions can transfer the 
Hanle rotation (of the plane of polarization) from the line core to the line wings before spontaneous 
de-excitation intervenes. In other words, in the presence of a small but significant clastic collision 
rate the Hanle effect does not vanish in the line wings. If the elastic collision rate is large then 
the collisions completely depolarize the scattered radiation throughout the line profile. This effect 
has been discussed in detail in HZ2 for the case ofaJ = 0->l->0 single scattering transition. 
However, these effects do not survive when the radiative transfer effects with angle-averaged PRD 
are explicitly taken into account. The collisional redistribution process is more effective in the 
case of angle-dependent PRD than in the case of angle-averaged PRD. Using the domain-based 
PRD theory of Bommicr [15] this effect was noticed even in the radiative transfer computations of 
Nagendra et al. [33] (see also Nagendra et al. [34]). It remains as effective in a pair of interfering 
doublet lines as in the case of a single line. However, in Sampoorna et al. [35] it was shown that 
the wing Hanle effect alone is insufficient to explain the observed wing signatures in the Q/I and 
U/I profiles of the Cai 4227 A line. 

As "f co h decreases to 0.5, which represents an equal mix of R 11 and R m , the values of Q/I in 
the wings of both the lines are significantly reduced (see dotted line in Figure 3). The collisional 
effects are now seen even in the core of the 5000 A line. This results in a decrease of the Q/I and 
U/I signals at the center of this line. The J-state interference signatures in Q/I are also modified. 




(50) 
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Figure 4: Emergent Stokes profiles at ^ = 0.047 computed for a slab of optical thickness T = 2 X 10 4 in the absence 
of a background continuum. The other model parameters are (a, e) = (10~ 3 , 10 -4 ). The magnetic field strength is 
set to zero. The coherence fraction is 7 co ^ = 0.1 (solid line), 7 co ^, = 0.5 (dotted line), y co h = 0.9 (dashed line), and 
"fcoh = 1 (dot-dashed line). 

When -fcoh is further reduced to 0.1, the effects of R m start to dominate over those of R n and also 
over the J-state interference effects (see dashed line in Figure 3). As a result the signatures of the 
J-state interference begin to fade away. The Q/I and U/I start to approach zero throughout the 
line profiles. As j co h is further decreased to 0.0001 (not shown in the figure), the collisional effects 
(through R m ) completely dominate the scattering process. This situation corresponds to a regime 
of extremely large line broadening. As a result the amplitude of / becomes much smaller compared 
to the other cases. Also, the Q/I, U/I, and V/I approach zero level throughout the line profiles. 

5.2. Polarized line profiles formed due to multiple scattering in an atmospheric slab 

In this section we present the emergent Stokes profiles computed by solving the polarized radia- 
tive transfer equation for a two-term atom including the effects of J-state interference and elastic 
collisions. For this we consider an isothermal constant property atmospheric slab with a given 
optical thickness T. The slabs are assumed to be self-emitting. The atmospheric model parameters 
used for the computations are represented by (T, a, e), where a is the damping parameter and e is 
the thcrmalization parameter defined in Equation (41) and the paragraph that follows. 

The Planck function B is taken as unity. The Doppler width for both the lines are assumed to 
be the same and equal to 0.025 A. For more details on the structure of the atmospheric slabs and 
the model parametrization we refer to P2. 

5.2.1. The non-magnetic case 

Figure 4 shows the emergent Stokes profiles which include the effects of J-state interference, 
clastic collisions and radiative transfer computed for a model atmosphere with parameters T = 
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Figure 5: Same as Figure 4 but computed for an optical thickness T = 2 X 10 8 . 

2 x 10 4 ,a = 10~ 3 and e = 10~ 4 in the absence of a background continuum. In these profiles the 
magnetic field is set to zero. Different line types represent different values of the coherence fraction 
jcoh- A range of values of ^ co h £ [1,0] is considered. As seen from the Figure 4, a decrease in "f C oh 
results in a gradual decrease in Q/I in the line core as well as in the PRD peaks of the 5000 A line. 
The intensity profiles are also quite sensitive to the effect of elastic collisions. As -y co h goes from 1 
(pure R n case) to 0.1 (R m dominated case), the self-reversed emission lines change over to nearly 
true absorption lines (thick solid lines). In Q/I the effects of collisions are confined only to the 
line core and the near wing PRD peaks. Specifically, it is shown by Nagendra [36] that the elastic 
collisions depolarizes the line core, and Te significantly depolarizes the line wing polarization 
(See Faurobert-Scholl [37]). The same conclusions are valid in the two-term atom model also. The 
interference region between the two lines seems to be less sensitive to the effect of elastic collisions. 

However for larger optical depths, significant dependence on 7 co ^ is exhibited in the wavelength 
region between the two lines. This can be seen in Figure 5 which shows the effect of elastic collisions 
in an optically thick atmospheric slab in the absence of a magnetic field. The model parameters arc 
the same as in Figure 4 but with T = 2 x 10 8 . As j co h decreases, the collisions take over the line 
formation process. When -y co h = 0.1 (thick solid line), deep absorption lines are formed in / with 
broad wings. The Q/I at the center of the 5000 A line becomes very small like in Figure 4. The zero 
crossing point at 5000.3 A remains the same for all the values of 7 co /i- In general a depolarization 
in Q/I is seen throughout the line profile because the radiative transfer effect is significant at all 
the frequencies. As expected, the Q/I reaches zero very far in the wings of both the lines after 
exhibiting a wing maximum nearly 10 A away from their line centers. The difference in behavior in 
the line core as well as in the line wings of the Q/I profiles formed under R n dominated (dot-dashed 
line) and R m dominated (thick solid line) conditions are better seen for the T = 2x 10 8 case when 
compared to the T = 2 x 10 4 case. 
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5.2.2. The magnetic case 

Figure 6 shows a comparison between the emergent Stokes profiles computed with (dashed 
line) and without (solid line) the presence of a weak magnetic field including the effects of elastic 
collisions. The magnetic profiles are computed for a field strength of w H = 0.004 with -f coh = 0.9. 
The model parameters are T = 2 x 10 4 , e = 10 -4 , a = 10~ 3 in the absence of a background 
continuum. An external weak magnetic field (through the Hanle effect) affects the multiply scattered 
Stokes profiles in a way similar to the singly scattered Stokes profiles. The Hanle effect causes a 
depolarization in Q/I at the center of the 5000 A line and also generates a U/I signal at this line. 
We recall that these effects are not seen at the 5001 A line since its polarizability factor Wi — 0. Like 
in the case of single scattered profiles, the magnetic field effects are confined only to the line core 
and the J-state interference signatures remain unaffected by the magnetic field. Also as discussed 
earlier, the wing Hanle effect in Q/I and U /I which were seen in the case of single scattered profiles 
in Figure 3 now disappear due to the radiative transfer effects. 

6. Conclusions 

In the present paper we have extended the theoretical framework for the J-state interference 
for type-II redistribution developed in PI, to include the effects of collisions (type-III redistribu- 
tion). The collisional PRD matrix is derived in the laboratory frame for a two-term atom with an 
unpolarized lower term and in the presence of magnetic fields of arbitrary strengths. However, the 
treatment is restricted to the linear Zeeman regime for which the Zeeman splitting is much smaller 
than the fine-structure splitting. The inelastic collisions coupling the upper term and the lower term 
and also the inelastic collisions coupling the fine structure states of the upper term are taken into 
account. However, the latter has been treated approximately. The approximation involves consid- 
ering only the depolarizing effects of the inelastic collisions but neglecting the polarization transfer 
rates between the fine structure states. The depolarization caused by the inelastic collisions has the 
same type of consequences as the depolarization by elastic collisions. Therefore we can merge both 
these effects into a common damping rate for the state Jb and appropriately redefine the branching 
ratios and the thermalization parameter for a two-term atom. A proper treatment of the inelastic 
collisions which cause polarization transfer requires formulating and solving the polarized statistical 
equilibrium equations. This is outside the scope of the present paper. The approximate treatment 
presented in this paper leads to slightly larger values of polarization in the line core as the inelastic 
collisions are not handled exactly. A treatment involving the statistical equilibrium equation would 
yield correct values of linear polarization. However in the line wings the formulation presented here 
becomes accurate enough and would give the same result as a full treatment in terms of statistical 
equilibrium equation, including PRD mechanism. 

The collisional frequency shift is inherently built into the redistribution matrix through the 
type-III redistribution function and the branching ratios. We discuss in detail the procedure of 
assigning the correct multipolar index K to the collisional branching ratio and depolarizing elastic 
collision rate D^ K ) . This procedure requires a detailed understanding of the role played by the 
multipolar index K for both the atom and the radiation field. We show how it becomes necessary 
to introduce a quantum generalized profile function for the case of a two-term atom in order to 
assign appropriate index K to the branching ratios and to D^ K \ In general is defined for 

each of the fine structure components (by making it depend on the quantum number Jb). However 
in the present paper we assume it to be independent of the J-quantum number. 
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Wavelength (A) 

Figure 6: Emergent Stokes profiles at /i = 0.047 computed for a magnetic field strength of »h = (solid line) and 
vji = 0.004 (dashed line) with a coherence fraction 7 co h = 0.9 in the absence of a background continuum. The model 
parameters are (T, a, e) = (2 X 10 4 , lO" 3 , 10~ 4 ). 
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Examples of the Stokes profiles resulting from single 90° scattering are illustrated for different 
values of the coherence fraction j co h- The profiles look similar to the ones presented in PI, which 
were computed using collisionless redistribution (the case of pure R n ), except for the depolarization 
in the wings of the Q/I profiles, and non-zero polarization in the wings of the U/I profiles. This 
interesting feature, which we refer to as the wing Hanle effect, is discussed. 

The effects of collisions are discussed by incorporating the newly derived collisional redistribution 
matrix in the polarized radiative transfer equation, in the simpler case of isothermal slab models. 
The technique of incorporating the Hanle redistribution matrix with the J-state interference and 
collisions, into the polarized radiative transfer equation for a two-term atom is presented. It is 
shown that the effects of elastic collisions in a two-term atom are similar to those of the two-level 
atom case. The redistribution matrices derived here have been used in the interpretation of the 
quantum interference signatures seen in the limb observations of the Cri triplet in Smitha et al. 
[10]. For simplicity the inelastic collisions between the fine structure states were neglected in that 
realistic modeling effort. 

With the present work we have further extended the theoretical tools that are needed for mod- 
elling the various spectral structures arising due to the transitions between fine structure states 
of an atom that have been observed in the Second Solar Spectrum so that they can be used to 
diagnose magnetic fields in regimes not accessible to the Zeeman effect. 
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Appendix A. Magnetic redistribution functions for type-Ill redistribution 



In this appendix we present the expressions for the magnetic redistribution functions of the type 
HH, HF, FH and FF appearing in Equations (16) and (17). They are defined as follows 
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and 



1 r+°° 

R^^b' f ( X b'f, x ba> ®) = 9 ■ rs I dlle 



«) 



a b + ( V 'b a ~ U ) 2 

sme' sine -—)■ (A ' 4) 

In the above equations H(a,x) and F(a,x) are the Voigt and Faraday- Voigt functions (see Equa- 
tion (18) of PI for their definition). is the scattering angle (the angle between the incident and 
scattered rays; see Figure 1 of PI). The dimensionless quantities appearing in Equations (A.l) to 
(A. 4) are given by 

VOba -v 76 +7c . , , v L 

Xba = — : ; a b = — — ; v bmllm = x ba + (g b fi b - g a ^a)-: , (A. 5) 

A^d 47tA^d Aj/d 

where v$ ba is the line center frequency corresponding to a J a — > J b transition in the absence of 
magnetic fields, a b is the damping parameter of the excited state 6, and Ai/q is the Doppler width. 
In the limit of a two- level atom (obtained by setting J b = J' b and J a = J/) and in the absence of a 
magnetic field the jj m > HH and R u ' u (defined in PI) reduce to the Hummer's R IU and R u functions 
respectively (see also HZ1). 

The angle-averaged analogues of Equations (A.1)-(A.4) are obtained through 

rf£™,b' m f m ( Xb 'f> X 'ba) = \J Rl b^,b'mfm( Xb 'f' X 'ba^ @ ) Sm0d0, (A.6) 

(see Equations (103) and (104) of Bommier [15] and Equations (30) and (31) of Sampoorna et 
al. [38]), where X and Y stand for H and/or F. A similar expression can be used for computing 
angle-averaged analogues of the type-II functions. 
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